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Abstract 

The extremely dismal prognosis of pancreatic cancer (PC) is attributed, at least in part, to lack of early diagnosis. Therefore, 
identifying differentially expressed genes in multiple steps of tumorigenesis of PC is of great interest. In the present study, a 
7,12-dimethylbenzanthraene (DMBA)-induced PC model was established in male Sprague-Dawley rats. The gene expression 
profile was screened using an oligonucleotide microarray, followed by real-time quantitative polymerase chain reaction 
(qRT-PCR) and immunohistochemical staining validation. A total of 661 differentially expressed genes were identified in 
stages of pancreatic carcinogenesis. According to GO classification, these genes were involved in multiple molecular 
pathways. Using two-way hierarchical clustering analysis, normal pancreas, acute and chronic pancreatitis, PanIN, early and 
advanced pancreatic cancer were completely discriminated. Furthermore, 11 upregulated and 142 downregulated genes 
(probes) were found by Mann-Kendall trend Monotone test, indicating homologous genes of rat and human. The qRT-PCR 
and immunohistochemistry analysis of CXCR7 and UBe2c, two of the identified genes, confirmed the microarray results. In 
human PC cell lines, knockdown of CXCR7 resulted in decreased migration and invasion. Collectively, our data identified 
several promising markers and therapeutic targets of PC based on a comprehensive screening and systemic validation. 
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Introduction 

Pancreatic cancer (PC) is a human solid malignant tumor with 
very poor prognosis [1]. Several clinical and pathological factors 
for PC have been identified, including T stage, lymph node/ 
distant metastasis, carbohydrate antigen (CA) 19-9 level and 
perineural/intraneural invasion; however, the factors affecting the 
prognosis of PC remain to be clarified [2,3,4,5]. In addition, 
molecular events involved in the pathogenesis and progression of 
PC, such as Kras mutation, have been discovered as hot spots [1]. 
However, more research is required identify the cellular factors 
that affect prognosis. 

One of the major drawbacks of the previous studies in 
human samples and cell lines is the difficulty to collect 
specimens in all stages of tumorigenesis. Therefore, animal 
models present unique advantages in this aspect. Chemical 
inducers of PC include N-nitrosobis(2-oxopropyl)amine, azase- 
rine, and 7,12-dimethylbenzanthracene(DMBA) [6,7,8,9,10]. 
These models can induce the entire range of carcinogenesis, 
through advanced stage and metastasis. DMBA has been widely 
used in the establishment of rat PC models 
[11,12,13,14,15,16,17]. In our previous studies, we found that 
acinar cells can transdifferentiate to ductal cells in the 
tumorigenesis process in the DMBA-induced PC rat model 
[18]. Nowadays, the notion and importance of acinar to ductal 



metaplasia (ADM) have been gradually accepted [19]. Previous 
studies have examined various features of these PC models, 
including histological/histochemical features [11,12,14], high- 
fat/high-protein diet as a promoter of carcinogenesis [13], 
glucose metabolism [15], and alterations in various proteins 
[16]. One recent investigation performed proteomic analysis in 
the PC rat model [1 7] . Thus far, the screening of differentially 
expressed genes in this model has not been reported. 

In the present study, our aim was to screen differentially 
expressed genes in PC using the DMBA-induced PC rat model. 

Materials and Methods 

Animals 

Adult Sprague-Dawley (SD) rats were provided by the 
Experimental Animal Center, Peking Union Medical College 
Hospital, Beijing, China. Rats were housed under standard 
conditions, including a pathogen-free environment and free access 
to food and drinking water. Twenty-four-hour urine samples were 
collected with metabolic cages in which only water but not food 
was provided. The Institutional Animal Care and Use Committee 
at Peking Union Medical College Hospital specifically approved 
this study and the use of rats. All efforts were made to minimize 
suffering. 
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Table 1. 


Mumbers of rats in different 


groups. 




Groups 


Time 


Number 


Experimenta 


7 days 


5 




2 weeks 


15 




1 month 


15 




3 months 


15 


Control 


7 days 


5 




2 weeks 


5 




1 month 


5 




3 months 


5 


Sham 


1 month 


5 
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Establishment of DMBA-induced PC Model in Rats 

The DMBA-induced PC model was established in male SD rats 
according to our previous method [14]. A total of 75 rats were 
divided into experimental, control and sham groups. DMBA or 
NaCl crystals (5 mg) were used for rats in the experimental and 
control groups, respectively, whereas no agent was applied in the 
sham group. In experimental and control groups, rats were 
sacrificed at 7 days, 2 weeks, 1 month and 3 months after 
implantation. Five rats in the sham group were sacrificed at 1 
month. Grouping of all the rats is shown in Table 1. 

RNA Extraction and Microarray Analysis 

Total RNA was extracted from -80°C frozen pancreatic tissue 
samples collected at all time points using the RNeasy mini kit 
(Qiagen, Germany) according to the manufacturer's instructions. 
Total RNA sample concentration and purity were examined and 
estimated by optical density measurements at 260/280 nm using a 
NanoDrop Spectrophotometer and agarose gel electrophoresis. 
The detailed measurement data and quality evaluation is shown in 
Table 2. The Superscript II reverse transcription kit (Invitrogen, 
USA), Genechip kit (Affymetrix, USA) and One-Cycle Target 
Labeling and Control Reagents (Affymetrix) were also used in the 
analyses. Hybridization was performed using the Genechip Rat 
Expression Set 230 (Affymetrix), according to the manufacturer's 
instructions. Gene expression analyses were performed using the 
Affymetrix Rat Genome 230 A array (Affymetrix, Santa Clara, 
CA, USA). Sequences used in the design of the array were selected 
from GenBank, dbEST and RefSeq. Gene chip 230 A contains a 
total of 15166 probe sets representing approximately 4700 known 
rat genes and 10460 unannotated expressed sequence tags. The 
accession number in GEO is GPL1355. 

Quantitative Reverse-transcription Polymerase Chain 
Reaction (qRT-PCR) 

Total RNA was extracted using TRIzol (Invitrogen). Reverse 
transcription was performed using the Superscript II reverse 
transcription kit (Invitrogen, USA). The TaqTM R-PCR SYBR 
Green I kit (Takara, Tokyo, Japan) was used for PCR. 
Amplification steps were as follows: 95°C for 30 s (pre-denatur- 
ation), and 40 cycles of 95°C for 5 s and 61°C for 31s. Beta-actin 
was designated as the internal control. All the experiments were 
repeated three times. Primers are listed in Table 3. 



Immunohistochemistry and Staining Evaluation 

Antibodies against CXCR7 and UBe2c were purchased form 
R&D and Abnova, respectively, and the PowerVisionTM two-step 
staining kit (PV-9000) was from Beijing Zhongshan Biotech Co., 
China. Briefly, sections were mounted, deparaffinized by xylene 
and rehydrated by ethanol. Antigen retrieval was performed by 
microwaving samples for 3 min. After blocking endogenous 
peroxidase, slides were then incubated overnight at 4°C with the 
primary antibody at a dilution of 1:200. Following washing in 
phosphate buffered saline (PBS), slides were incubated with 
horseradish peroxidase (HRP)-labeled secondary antibody for 
30 min at room temperature. Diaminobenzidine was used as a 
chromogen. Finally, slides were counterstained with hematoxylin. 
Brown coloration in the cytoplasm or nucleus was defined as a 
positive signal. The ratio percentage of staining positive stained 
cells more than 50% was classified as a the strong positive signal. 

Cell Culture and Knockdown of CXCR7 

Two human pancreatic cancer cell lines (PANC-1 and SU86.86) 
were kind gifts from Professor Helmut Freiss, Heidelberg 
University, Germany [20,21]. Cells were cultured in RPMI 1640 
medium and DMEM medium, with 10% fetal bovine serum 
(Invitrogen) at 37°C in a humidified atmosphere of 5% C0 2 . Cells 
seeded in 6-well plates were transfected with CXCR7 or control 
(scramble) small interfering RNA both at a concentration of 
40 nM using Lipofectamine RNAiMax (Invitrogen), according to 
the manufacturer's instructions. The CXCR7 RNAi was pur- 
chased from Invitrogen and the sequence is as follows: AGAU- 
GUAGCAGUGCGUGUCAUAGCC. 

Western Blotting 

Cells were washed with PBS and collected by scraping into 
individual tubes. Proteins were extracted according to protein 
extraction protocols, and protein concentrations were determined 
using the Pierce BCA protein assay kit (Thermo Scientific, 
Meridian Rd, Rockford). Protein extracts (80 ugVlane) were 
electrophoresed on 10% polyacrylamide gels (SDS-PAGE) fol- 
lowed by transfer to PVDF membranes and blocking with 5% 
non-fat dry milk for 2 h. Membranes were incubated with primary 
antibody against CXCR7 (Abgent, San Diego, CA) overnight at 
4°C. Secondary antibody (anti-rabbit IgG) was incubated for 1 h 
at 37°C. Blots were washed three times with PBS, exposed to 
chemiluminescence reagents (Merck Millipore, Darmstadt, Ger- 
many), and exposed to photographic film. All experiments were 
performed in triplicate. 

Cell Migration and Invasion Assays 

Transwell inserts (8.0 |im pore size) (CORNING, Chelmsford 
St) were used for migration and invasion assays. The lower 
chamber was filled with 500 |J.L of RPMI 1640 medium or 
DMEM medium with 10% FBS. PANC-1 and SU86.86 cells were 
resuspended in the migration medium (serum-free RPMI 1640 or 
DMEM) and 4xl0 5 PANC-1 or SU86.86 cells in 200 uL of 
medium were added to the upper chamber after being washed 
twice with PBS. After 24 h of incubation at 37°C, the cells on the 
upper surface of the membrane were gently and carefully scraped 
out using cotton tips. The migrant cells that were adherent to the 
lower surface of the membrane were fixed in 10% formalin at 
room temperature for 30 min, and then stained with H&E 
(hematoxylin and eosin stain). The results of the cell migration 
assay were evaluated by counting the number of cells on the lower 
surface of the membrane under an inverted microscope in five 
different fields at a magnification of x200. 
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Table 2. Total RNA measurement data and quality evaluation. 





Sample Serial 
Nunber 


RNA Concentration (ng/fi D 


A 2 6o/A 2 8o 


RNA Quantity (fig) 


Electrophoresis 


Final Evaluation 


Norm-panl 


2.448 


2.03 


489.5 


Unqualified 


Unqualified 


D4001 


0.905 


2.03 


181.1 


Qualified 


Qualified 


D4002 


1.776 


2.03 


355.1 


Unqualified 


Unqualified 


D4003 


1.375 


2.07 


275.0 


Qualified 


Qualified 


D4004 


1.910 


2.06 


381.9 


Unqualified 


Unqualified 


DB7001 


0.917 


2.04 


183.3 


Qualified 


Qualified 


DB7002 


2.882 


2.02 


576.3 


Unqualified 


Unqualified 


DB7003 


3.511 


1.96 


702.1 


Qualified 


Qualified 


DB7004 


1.872 


2.05 


374.4 


Qualified 


Qualified 


D7001 


1.849 


2.05 


369.9 


Qualified 


Qualified 


D7002 


2.344 


2.04 


468.8 


Unqualified 


Unqualified 


D7003 


1.932 


2.05 


386.4 


Unqualified 


Unqualified 


D7004 


0.542 


1.98 


108.5 


Qualified 


Qualified 


D7005 


2.168 


2.03 


433.5 


Qualified 


Qualified 


WB2001 


2.229 


2.03 


445.7 


Unqualified 


Unqualified 


WB2002 


3.356 


1.99 


671.3 


Unqualified 


Unqualified 


WB2004 


1.532 


2.05 


306.3 


Qualified 


Qualified 


WB2005 


0.437 


2.01 


87.4 


Qualified 


Qualified 


W2001 


1.175 


2.03 


235.1 


Qualified 


Qualified 


W2004 


2.604 


2.00 


520.8 


Qualified 


Qualified 


W2005 


1.539 


2.03 


307.9 


Qualified 


Qualified 


W2006 


2.451 


2.02 


490.2 


Qualified 


Qualified 


W2007 


1.829 


2.01 


365.7 


Qualified 


Qualified 


W2009 


2.660 


2.00 


531.9 


Qualified 


Qualified 


W2010 


1.425 


2.03 


285.0 


Qualified 


Qualified 


W2011 


0.442 


2.00 


44.2 


Qualified 


Qualified 


W2012 


1.949 


2.04 


389.8 


Qualified 


Qualified 


W2013 


1.413 


2.02 


282.6 


Qualified 


Qualified 


W2014 


2.270 


2.02 


453.9 


Qualified 


Qualified 


MB1001 


3.368 


1.95 


673.7 


Qualified 


Qualified 


MB1002 


1.449 


2.00 


289.7 


Qualified 


Qualified 


MB1003 


2.696 


1.99 


539.2 


Qualified 


Qualified 


MB1004 


2.934 


1.98 


586.7 


Unqualified 


Unqualified 


Ml 005 


3.510 


1.94 


701.9 


Qualified 


Qualified 


Ml 01 0-2 


0.970 


1.99 


194.0 


Qualified 


Qualified 


M1012 


1.387 


2.02 


277.5 


Qualified 


Qualified 


M1013 


1.727 


2.04 


345.4 


Qualified 


Qualified 


M1014 


1.983 


2.01 


396.6 


Qualified 


Qualified 


M1015 


2.701 


1.99 


540.2 


Qualified 


Qualified 


M3001 


3.039 


1.98 


607.9 


Qualified 


Qualified 


M3002 


1.624 


2.02 


324.9 


Qualified 


Qualified 


M3004 


0.891 


1.98 


89.1 


Qualified 


Qualified 


M3005 


1.437 


2.02 


287.4 


Qualified 


Qualified 


M3006 


0.432 


2.02 


21.6 


Qualified 


Qualified 


M3008 


0.392 


2.02 


19.6 


Qualified 


Qualified 


M3009 


1.985 


2.02 


397.0 


Qualified 


Qualified 


M3014 


1.466 


1.99 


293.2 


Qualified 


Qualified 
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Table 3. Primers used in this study. 





Gene 


Primers 


length (bp) 


P-actin 


5'- TCTGTGTGGATTGGTGGCTCT -3' 
(forward) 


300 




5'- AGAAGCATTTGCGGTGCAC -3' 






(reverse) 




CXCR7 


5'- GAGCATC I I LI ICCTCGCATG -3' 
(forward) 


300 




5'- CAACACGGCGTACCATCTTCT -3' 






(reverse) 




ATP6v1g2 


5'- GAATCTCTGCGTCCCATGACA -3' 
(forward) 


250 




5'- CCCCATCACCACGACATACAT -3' 






(reverse) 




UBe2c 


5'- TTCAAAGCAGGTCTCCAACCA -3' 
(forward) 


400 




5'- CAGCTCAGAAACCACGGAGAA -3' 






(reverse) 
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For cell invasion assays, the undersurface of the membrane was 
coated with ECM gel (Sigma-Aldrich, Saint Louis, MO, USA) in 
PBS for 5 h at 37°C. PANC-1 or SU86.86 cells (8xl0 5 ceUs) in 
serum-free RPMI 1640 medium or DMEM medium were seeded 
onto ECM gel-coated upper chambers, according to the manu- 
facturer's instructions. Cells that had crossed the ECM and passed 
through the pores of the filter were counted in five fields at x200. 
The results are representative of three different experiments. 

Statistical Analyses 

For microarray results, background correction and normaliza- 
tion were performed. One-way analysis of variance (ANOVA) and 
Student-Newman-Keuls (SNK) tests were used to detect statisti- 
cally differentially expressed genes between different groups. In 
bioinformatic analysis, differentially expressed genes were subject- 
ed to GO analysis hierarchical clustering, principal component 
analysis, SOM and Mann-Kendall test. The Student-t and Chi- 
square tests were applied to show differences of measurement and 
categorical data. The statistical software package SPSS 13.0 (SPSS 
Inc., Chicago, 111) was used. A statistically significant P value was 
defined as P<0.05. 




Figure 1. Pancreatic tumor formation and metastasis in the DMBA-induced PC rat model. (A) 

Pancreatic tumor of 3-month group. (C) Small intestine metastasis nodules of pancreatic tumor. (D) Liver 
doi:1 0.1 371 /journal.pone.008291 0.gOOl 



Pancreatic tumor of 1 
metastasis nodules of 



-month group. (B) 
pancreatic tumor. 
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Figure 2. Histological changes in different stages of the DMBA- 
induced PC rat model. (A) Normal pancreas (original magnification 
x60). (B) Acute pancreatitis (original magnification x150). (C) Chronic 
pancreatitis (original magnification x150). (D) PanIN 1a (original 
magnification x60). (E) PanIN 2 (original magnification x150). (F) PanIN 
3 (original magnification x60). (G) Well-differentiated PC (original 
magnification x60). (H) Poorly-differentiated PC (original magnification 
x60). 

doi:1 0.1 371 /journal.pone.008291 0.g002 

Results 

Successful Establishment of DMBA-induced PC Model in 
Rats 

Acute pancreatitis was observed in rats (7 days after implanta- 
tion) in both the experimental and control groups. Pancreatic 
intraepithelial neoplasia (PanIN) was present in 6 out of 15 rats 
(40%) in the experimental group 2 weeks after implantation, 
whereas no PC was found in rats in the control group. One month 
and 3 months after implantation, 12 out of 15 (80%) and 14 out of 
14 (100%) rats in the experimental group, respectively, developed 
PC (Fig. 1A, IB). Four rats in the 3-month group developed small 
intestine metastasis or liver metastasis (Fig. 1C, ID). However, PC 
was absent in rats in the control groups. The PC formation rate in 
the experimental group was higher than that in the control group 
(P<0.05). One rat died before the endpoint (3 months), causing a 
mortality rate of 6.7% (1/15). 

Histological changes were evaluated using HE staining and light 
microscopy. Normal pancreas (Fig. 2A), acute pancreatitis 
(Fig. 2B), chronic pancreatitis (Fig. 2C), PanIN 2a (Fig. 2D), 



PanIN 2 (Fig. 2E), PanIN 3 (Fig. 2F), well-differentiated PC 
(Fig. 2G) and poorly-differentiated PC (Fig. 2H) were observed. 

Differentially Expressed Genes in Rat DMBA-induced PC 
Model 

A total of 661 differentially expressed genes were detected in 
normal pancreas from the control group and the experimental 
groups (Table 4). According to GO classification, the genes were 
involved in various functional categories and multiple biological 
processes, including gas transport, carbohydrate metabolic pro- 
cess, antigen processing and presentation (Fig. 3A). By two-way 
hierarchical clustering analysis, normal pancreas, acute and 
chronic pancreatitis, PanIN and early and advanced pancreatic 
cancer samples were completely discriminated (Fig. 3B). There 
were 11 upregulated genes (probe) and 142 downregulated genes 
(probe) by Mann-Kendall trend Monotone test (P<0.05) (Fig. 3C). 
The top 20 upregulated and downregulated genes in the DMBA 3 
month group whose expression levels were changed by two-fold 
compared with the control group (P<0.01) are shown in Table 5. 
In addition, homologous genes in rat and human at different stages 
of the rat model were successfully screened (Table 6, Fig. 3D). 

Validation of Differentially Expressed Genes 

We selected three of the differentially expressed genes, CXCR7, 
ATP6vlg2 and UBe2c, for further examination. First, we 
examined their expression profiles by qRT-PCR. The ACt values 
of CXCR7 gradually decreased, along with prolonged DMBA 
treatment (P<0.05) (Table 5), indicating upregulated expression. 
There were no significant differences for ATP6vlg2 and UBe2c 
(P>0.05) (Table 7). Using immunohistochemistry, through a 
pathologist's diagnosis and confirmation, significantly enhanced 
CXCR7 and UBe2c expression were detected, along with the 
progression of PC in rats (P<0.05) (Table 8, Figs. 4). 

Decreased Migration and Invasion after Transfection of 
CXCR7 siRNA in PANC-1 and SU86.86 Cell Lines 

Expression of CXCR7 was successfully downregulated in the 
two human PC cell lines, PANC-1 and Su86.86 (Fig. 5A), after 
transfection with CXCR7 siRNA, compared to transfection with 
control siRNA. Significant differences in cell migration and 
invasion were observed in CXCR7 siRNA transfected cells 
compared to control cells (P<0.05) (Figs. 5B and 5C). 

Discussion 

Various alterations in molecular pathways have been known to 
play an important role in the initiation and progression of PC [1]. 
However, one of the complications in studies based on human 
samples is that features of each disease phase cannot be easily 
displayed. In this aspect, animal models carry a unique advantage. 
In the field of PC research, major models have included DMBA- 
induced and transgenic models [11,12,13,14,15,16,17,22]. In 
studies of the DMBA-induced rat model, different DMBA doses 
and rates of PC development have been reported 
[11,12,13,14,15,16,17]. The current study, which was based on 
our previously reported method [14], achieves satisfactory PC 
induction results. Therefore, these findings indicate successful 
establishment of the model. 

Microarrays have been widely applied in the screening of 
differentially expressed genes in PC [23,24,25,26]. However, this 
technique has not been used in the DMBA-induced PC model. 
The present study examined the gene expression profile in this PC 
rat model. Affymetrix microarray analysis identified a total of 66 1 
genes that were differentially expressed in normal pancreas from 
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Figure 3. Microarray results. (A) GO classification showing genes involved in different molecular function categories. (B) Two-way hierarchical 
clustering analysis. (C) P values in Mann-Kendall trend Monotone test. (D) Homologous genes in rat and human at different stages of the rat model. 
doi:10.1371/journal.pone.0082910.g003 



the control group and the experimental groups. These genes 
included critical facilitators of important functions during pancre- 
atic tumorigenesis [27], such as Kras, CDKN2A, SMAD4, and 
TP53. These genes were upregulated or downregulated during 
pancreatic tumor formation duration. GO classification, two-way 
hierarchical clustering analysis and Mann-Kendall trend Mono- 
tone test further identified genes, including nm23, cyclin Bl and 
FGFR3. A previous study showed that nm23 was associated with 
metastatic potential of human PC cell lines [28]. In the current 
study, we showed that nm23 expression was significantly different 

Table 4. Differentially expressed genes in different groups of 



rats (SNK test). 






Normal 


DMBA-7d 


DMBA-2W 


DMBA- 
lm 


DMBA- 
3m 


Normal 




1765* 


2815* 


4826* 


2435* 


DMBA-7d 


1765* 




413* 


1002* 


566* 


DMBA-2W 


2815* 


413* 




530* 


393* 


DMBA-lm 


4826* 


1002* 


530* 




947* 


DMBA-3m 


2435* 


566* 


393* 


947* 




*P<0.05. 

doi:1 0.1 371 /journal.pone.008291 0.t004 



between normal pancreas and the DMBA 1 month group, 
suggesting that the molecule might also be involved in the early 
stage of pathogenesis of PC. Little data is available concerning the 
roles of cyclin Bl and FGFR3 in PC [23,29]. Thus here, we 
provide further evidence that substantiates their involvement in 
this disease. Further investigations on the detailed mechanisms 
underlying their involvement in PC are required. 

Expression patterns of three differentially expressed genes, 
CXCR7, ATP6vlg2 and UBe2c, were evaluated by qRT-PCR 
and immunohistochemical staining. The results showed gradually 
upregulated CXCR7 expression at the mRNA level, and 
enhanced positive expression of CXCR7 and UBe2c, along with 
prolonged DMBA treatment. These findings mosdy followed the 
microarray data. CXCR7 was previously demonstrated to have an 
impact on growth, apoptosis, invasion /metastasis and prognosis in 
many cancers [30,31,32,33]. A recent study revealed that CXCR7 
promoted cell growth in PC [34]. However, its prognostic 
significance in PC remains controversial [35,36]. In the rat model 
and human cell lines used in our experiments, the association of 
CXCR7 and progression as well as invasive proclivity of PC was 
indicated. Therefore, further mechanistic studies for the molecular 
factors involved in PC worth are warranted. 

UBe2c was established to be associated with growth and 
apoptosis inhibition in cancer cells, thus playing a role in 
tumorigenesis and showing potential to be a biomarker of both 
diagnosis and prognosis [37]. However, UBe2c was not previously 
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Table 5. List of upregulated and downregulated genes in DMBA-3 month group whose expression levels were changed by more 
than 2-fold compared with the control group {p<0.01). 





Top 20 Genes 


Upregulated: 


CXCR7, Clcn3, Kras, cyclinBI, nm23, ATP6v1g2, SMAD4, Acf, Adam6, Cyb5b, FGFR3, Mndl, Kifcl, Car8, Ptpla, EGFR, Park2, c-Myc,Cadps2, 
Etv5 


Downregulated: 


Emb, Gabbrl, Zdhhc8, CDKN2a, CXCL10, TP53,Gbp2, Tnfsf9, Dnaselll, Fgr, Vcan, Lyz,Plod2, Stat5a, Smarca2, Cspg4, Rgs14, Ldb3, Bcl2a1, 
CCR5 
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Table 6. Homologous genes in rat and human at different stages of rats (SNK test). 
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*P<0.05. 
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Figure 4. Expression of CXCR7 and UBe2c in DMBA-induced PC rat model (3-month) and normal rat pancreas in control groups. 

(original magnification x200). (A),(C) Normal rat pancreas,HE staining. (B),(D) Rat PC : HE staining. (E),(G) Normal rat pancreas,CXCR7 
immunohistochemistry staining. (F),(H) Rat PC,Ube2c immunohistochemistry staining. 
doi:1 0.1 371 /journal.pone.008291 0.g004 
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Table 7. Validation of some differentially expressed genes by qRT-PCR (ACt value). 





Genes 


Normal 
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DMBA-2W 


DMBA-lm 
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P value 


UBe2c 
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8.26±0.88 
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11.45±0.46 


0.574 


CXCR7 


7.65 ±0.65 


7.45 ±0.1 6 


7.10±0.87 


6.56±0.22 


6.06±0.29 


0.046 
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Table 8. Validation of some differentially expressed genes by Immunohistochemistry. 
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Figure 5. Knockdown of CXCR7 resulted in decreased migration and invasion in human PC cell lines. (A) Knockdown of CXCR7 by siRNA 
transfection. (B) Migration after transfection of CXCR7 or control siRNAs. (C) Invasion after transfection of CXCR7 or control siRNAs. 
doi:1 0.1 371 /journal.pone.008291 0.g005 
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investigated in PC. In the present study, we found enhanced 
positive expression of UBe2c along with prolonged DMBA 
treatment, suggesting its potential role in progression of PC. 
Notably, UBe2c mRNA expression was not significandy altered, 
indicating that post-transcription regulation might be involved in 
changes of its expression. Detailed mechanistic data underlying 
this activation is required. 
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